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Abstract 

In this paper we present results of computations for the ground energy of 
weakly coupled double well potential in quantum mechanics. We give a numer- 
ical evidence for cancelation of imaginary contributions to energy coming from 
Borel resummation and multi-instanton terms. We also estimate several higher 
coefficients of the multi-instanton expansion which are not given in the literature. 

1 Introduction 

The question of the relation between the perturbative and non-perturbative con- 
tributions in quantum theories is a long standing one. The answer emerged over 
years of studies and is well known, e.g. for the anharmonic oscillator [UH]. The 
two contributions are in principle additive, however there is a sublte interplay be- 
tween them. Namely, the ambiguities in resummation of perturbative, asymptotic 
series are non-perturbative and cancel against ones of multi-instanton contribu- 
tions. In this paper we verify this claim by confronting theoretical predictions 
with, very high precision, numerical solutions of the problem. 

We consider a double well potential in the following parametrization: 

V{x) = \x 2 {l- y fgxf. (1) 

As it is well known [3], the perturbation series for the ground energy E(g) is 
not Borel summable. Still, one can perform a Borel sum for complex coupling 
constant g and analytically continue it to positive axis from lower or upper half 
of the complex plane. This freedom results in an ambiguity of imaginary part of 
the Borel sum. However, there is another ambiguity originating from instanton 
contributions to energy. As stated in [I] these imaginary terms must cancel. 

In this paper we give a direct numerical evidence of this fact. To this end, 
we find perturbative series of the ground state energy and construct its Borel 
transform. Then, we continue it analytically using Pade approximation and 
perform inverse Borel transform. It is done in the limit g = Keg + i0 + . Next, 



* zbigniew. ambrozinski@uj . cdu . pi 
twosiek@th.if.uj.edu.pl 



we demonstrate that imaginary part of the two-instanton molecule contribution 
derived by Bogomolny [2] cancels imaginary part of energy in Borel sum at leading 
order. We also confirm cancelation of higher order imaginary terms given by 
Jentschura and Zinn-Justin in [?] and find next few coefficients of their expansion. 

To compute energies numerically we use, essentially exact, cut Fock space 
approach [S] H3 [Hj . With this method we show that real part of the Borel sum 
plus two-instanton molecule contribution give very accurate approximation to 
the energy for small couplings. 

2 Borel resummation 

One can find perturbation series of the ground state energy E(g) = YlkLo e kg k 
up to high orders using the Rayleigh-Schrodinger perturbation theory [5]. We 
found £fe for k < 500. Asymptotic behavior of ej, is known [10] and yields 

e fc ^-fc!3 fe -. (2) 

7T 

The relative difference between asymptotic estimate and exact value of is 0.6% 
for k = 500 and decreases at rate estimated in [TJjj. The perturbative series is 
asymptotic and a resummation procedure is needed. To this end we use the Borel 
transform 

fc=0 

■Boo(t) is convergent for |t| < | and has a pole at t = |. Both follow from the 
asymptotic behavior of e&. Inverse Borel transform is given by the integral 

E Borel {g) = - dte' t/9 B(t) (4) 
9 Jo 

where B(i) is analytic continuation of the Borel transform Baa(t). Because the 
pole at t = | lies on the integration path, the perturbative series is called not 
Borel summable [3]. Still, the integral can be calculated for t — Ret ± i0 + . 
One can show that it is equivalent to take t real but g — Re<? ± i0 + . Changing 
between positive and negative imaginary part of g alters sign of imaginary part 
of the integral Q. In particular, contribution of the leading singularity at t = | 
is 

E Borel (g + i0 + )-E Borel (g-i0 + ) = -^ e - 1/3g . (5) 

As an approximation of the analytic continuation B(t) we took the diagonal 
/ P ac ^ approximant P(t) of Bxit). We shall now analyze poles of P(t). 
They are presented in Fig. [l] Poles on the real axis condense with growing K 
and form a cut on the interval [|,oo). They reflect a cut of B(t). Poles with 
nonzero imaginary part move to infinity as K increases. We infer that they lie in 
the region where the Pade approximant is no longer reliable. It is known [11] that 
the approximant is weakly convergent near poles of the approximated function. 
Therefore, we changed the integration contour from t = Ret + i0 + to t = |t|e l7r / 4 . 
Secondly, we cut the integral at Ret = 1.4 so that we did not come close to poles 
of P(t). Error coming from both, change of integration contour and cutting the 
integral is of order ^e~ 1A ^ 9 . It is much smaller than the ambiguity of Borel 
sum already for g = 0.1 and therefore we shall neglect it. Major error for small 
coupling constant g is an effect of finite K. It is a nontrivial task to estimate it 
a priori. One has to try different K's and check if one can reach such values that 
the energy is independent of K. 



2 



Re Re Re 

(a) K = 200 (b) K = 350 (c) K = 500 

Figure 1: Poles of Pade approximant for different orders of approximation K. 



3 Instanton contributions 

The perturbation theory yields the same series expansion of the ground energy 
in both minima of the potential V(x), x = and x = l/^g- Thus, the ground 
energy is degenerate at the level of perturbation theory. 

One can calculate splitting of the ground energy Eq and the first excited 
energy Ei using semiclassical methods in the dilute instanton gas approximation. 
The difference is E^ — Eq = ^=e _1/ ^ 6s . The leading correction to dilute gas 
approximation due to interactions between instantons was later calculated by 
Bogomolny [2J: 

E%\g) = ^e- 1 ' 33 ( j + ln (-?)), N = 0,1, (6) 



^9 V V Sjj 

where 7 w 0.577 is the Euler constant. One has to understand this formula as a 
continuation from negative g through upper or lower half of the complex plane, 
i.e. taking g + i0 + or g — i0 + limit. These limits give different results: 

E^\g + tQ + )-E^(g-iQ + )^ 2 je- 1 ^, N = 0,1. (7) 

Note that this is exactly opposite to the leading order of ambiguity of Borel sum 

Full formula for instanton and perturbative contributions to energy was given 
by Zinn- Justin Q] [4] : 

00 00 , „-l/6g\ nn - 1 / / o\\i 00 

^( ff ) = E^ fc + E(-(-^^r) EM" ^ (8) 

fc = n=l ^ V » / ; = V V »// fc=0 

for = 0,1. Some coefficients were given in [3]: 

£200 = 7, £210 = 1, 

23 53 53 
£20i = -y-y7, 6211 = —, (9) 

13 1277 1277 

£202 = -= ^~7i £212 
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Cancelation of ambiguities (|5| and Q renders formula (JsJ> unique at least at 
order g~ 1 e~ 1 ^ 3a . We will use numerical analysis to see that the series is unique 
also at higher orders in g, i.e. n = 2 and k > 0. We will also show that real part 
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Figure 2: Convergence of eigenvalues of the matrix Hm with growing M for g — 0.002. 
Behavior of eigenenergies changes qualitatively around the top of the barrier. 

of energy improves when one adds two-instanton terms to the Borel energy 
This can be done only if we eliminate much larger contributions from independent 
instanton (n — 1). Note that n = 1 terms are exactly opposite for the ground 
and first excited energy. Therefore, we will be interested in their mean values 
E=±(E +E 1 ) 

4 Cut Fock space method 

An alternative technique of computing the lowest energies is the cut Fock space 
method. Let us denote by \n) the Fock basis, which is set of energy states of 
the harmonic oscillator with minimum at x — tt^. Then the matrix Hm = 

((n||P 2 + V(X)\m)) nrn with n,m < M is an approximation of the Hamiltonian. 
M is called the cut-off. We expect that the lowest eigenvalues of Hm approximate 
energies of the system. It was very well confirmed in many cases [5l El E] - For 
our system convergence of energies with increasing M is presented in Fig. [2] 

For small coupling constant g there are several eigenenergies smaller than 
height of the barrier, e.g. there are 38 for g — 0.002 with V(l/2y/g) « 15.6. 
Classically the potential barrier is impenetrable for states with such energy. In 
the classical limit, there is a pair of degenerate states, one localized in left and 
one in right minimum of the potential. Due to quantum tunneling these states 
mix into symmetric and antisymmetric combinations. The symmetric state of 
each pair has slightly lower energy than the antisymmetric one. 

Let us now analyze convergence of the energies with growing cut-off. Char- 
acteristic feature of this method is, as in any variational method, that energies 
are approximated from above. Therefore, they fall down as the cut-off increases 
and approximation improves. For small cut-offs energies are higher than the po- 



4 



tential barrier and splitting between them is large. Only when energies become 
smaller than the potential barrier they join into pairs. This can be clearly seen 
in Fig. [| 

Another feature is that energies fall linearly for small cut-offs. This is because 
the space basis does not yet explore minima of the potential. For cut-offs high 
enough the convergence becomes approximately exponential. One can find a 
more detailed analysis of this system in |12| . 

The cut Fock space method yields no approximations apart from precision 
of computations and finite cut-off effects. Thus, it is essentially exact. For our 
purposes we consider the value 

E Fo ck{g) = \{E + Ei) (10) 

where _Eo,i are the two lowest eigenvalues of Hm- For the smallest considered 
value of coupling constant g = 0.00016 the needed cut-off was M — 5000 and 
precision 10 -980 . The cut Fock space method turns out to be far more efficient 
than computing energies through Borel resummation procedure. Still higher 
cutoff and greater numerical precision can be applied resulting in more accurate 
results. 



5 Comparison 

We will show that ImE Borel (g) + Im_B (2) (g) = and ReE Borel (g) + Re£ (2) (g) = 
EFockig) up to leading order disregarded in E^ 2 \g). To this end we introduce 

^^^tM-llEW, (id 

y i=o v V a/ / fc=0 

K _ lmE Borel (g)+lmE^ K (g) 

Al W - I e -l/3 9 ' <- 12 ) 

3 

a k, v ReE Borel (g) + ReE^' K (g) - E Fock (g) 

Ani9) - ^-^H2/g) • (13) 

Denominators of (g) and A^(g) are leading terms of imaginary and real part 
of E^ 2 \g) respectively. According to formula |8j) we expect that A^ x (g) = 
0(g K+1 ). Plots of x (g) for K = 0, 1, 2 are presented in Figs 3ji Asymptotic 
behavior of A^ x (g) agrees with predictions. We also found approximations 
of coefficients Q from numerical data and confirmed coefficients £2ife in with 
precision 10 -20 and coefficients £20fe with precision fO -8 . 

It turned out that it is possible to determine several next coefficients t2ik in 
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Figure 3: Function Ag(g) for if = 0,1,2. One expects that Ag(g) = 0{g aK ) with 
= K + 1. Fitting a straight line to the log-log plot gives a K =o = 1-0242 ± 0.0007, 
a K= i = 2.05 ± 0.02, a K=2 = 3.039 ± 0.002. Fitted values of a K are slightly above 
expectance because of higher order corrections, which are not negligible for nonzero g. 
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Figure 4: Function A%(g) for K = 0,1,2. One expects that A%(g) = 0{g aK ) with 
ax = K + 1. Fitting a straight line to the log-log plot gives cck=o = 1-001 ± 0.0001, 
a K =i = 2.0075 ± 0.0008, a K=2 = 3.009 ± 0.0009. 
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expansion Q. Taking coefficients (JoJ) as given we found 
45941 336437 1n _ xo 

£203 = 7 ± 1.0 • 10 , 

144 1296 ' 

20772221 141158555 . „ ,„_ 6 

£204 = 7 ±2-10 , 

2592 31104 ' 

£ 205 = -205496.5847- ^£^7 ± 2 • 10" 3 , 
186624 

e 206 = 69 3 6 9 80.4 ± 4.8 

336437 21 
6213 = -^296" ±L3 ' 10 ' 



£214 



141158555 
31104 



17^42610737 13 

£215 = ± 5.9 • 10 , 

186624 

£216 = -2221191.7314262645 ± 4.8 • 10" 9 , 
£217 = -58524267.633067 ± 2.5 • 10~ 5 , 
£218 = -1695080020.213 ± 9.5 ■ 10" 2 , 
£219 = -53461315700 ± 1.6 • 10 3 , 
£2i ; io = -1823771270000 ± 4.8 • 10 5 . 

In Fig. [H] the plot of Ax(ff) with found coefficients is presented. One can 
observe that A% (g) decreases with growing K and the relation Af (g) = 0{g K+1 ) 
holds. Having already a few coefficients of the n = 2 expansion we found its 
Borel sum. It turns out that this procedure improves convergence by at least 
a factor of 10 2 for g £ (0.00016,0.009). From formula (JSJl it can be seen that 
n — 2 terms do not have to remove full ambiguity of the energy. Some of it may 
be removed by n = 4 contribution. Therefore, we cannot expect Ar_(g) to be 
smaller than ^e^ 1 ^ 39 ln 2 (|). This limitation is indeed seen in Fig. 5 Similar 

analysis was made for Ax (g) and results are shown in Fig. |6] A detailed analysis 
concerning cancelation of higher order ambiguities was presented in [13] for the 
cosine potential. 



6 Summary 

Using essentially exact numerical solutions of the anharmonic oscillator problem, 
we have verified the long-existing theoretical predictions that the ambiguities 
in resumming the perturbative, asymptotic series are nonperturbative, and in- 
deed cancel with the ones from two-instanton interactions. The remaining, well 
defined part, is constructed additively from resummed perturbative series and 
nonperturbative contributions and is in agreement with our numerical results for 
g < 0.05. For agreement in higher orders the story repeats itself on the level of 
multi-instanton interactions. 

Thanks to high precision of computations we were able to confirm values of 
coefficients of two-instanton correction given by Zinn- Justin and Jentschura. It 
was also possible to estimate a few more coefficients of the energy expansion 
which appear in imaginary part of energy. Computing energies for smaller values 
of the coupling constant g would require yet higher precision. For g = 10~ 5 the 
instanton correction is of order 10 -14471 so precision of computations would have 
to be 16 times higher than for g = 0.00016. 

The series E^ 2 '{g) appears to be asymptotic and not Borel summable. Con- 
siderable improvement of results is observed when on finds Borel sum of E^ 2 \g) 
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Figure 5: Dashed lines represent function A% (g) for K = 2, 4, 7, 10. Imaginary ambi- 
guity of energy decreases for growing K. Thick line is Aj (g) with instanton contribu- 
tion replaced by its Borel sum. Solid line is the leading term of the n = 4 contribution 
in expansion ([8]) with £430 = 1. Since it was neglected in our analysis, it is the lower 
bound for Atf(g). It is approximately saturated for g > 0.02. 
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Figure 6: Dashed lines represent function A^(g) for K = 0,2,4,6. Approximation 
of energy improves for growing K. Thick line is A^(g) with instanton contribution 
replaced by its Borel sum. Thin line is n — 4 contribution which is neglected in our 
analysis. 
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even with only few coefficients of given expansion. Having more terms it might 
be possible to extract (g) contribution to the ground energy. 

Acknowledgements 

We thank Mithat Unsal, Michael Teper and Yannick Meurice for fruitful discus- 
sions. This work was supported by Foundation for Polish Science MPD Pro- 
gramme co-financed by the European Regional Development Fund, agreement 
no. MPD/2009/6. 

References 

[1] J. Zinn- Justin, Nucl.Phys. B192 (1981) 125-140. 

[2] E.B. Bogomolny, Phys.Lett. B91 (1980) 431-435. 

[3] E. Brezin, J.-C. Le Guillou, and J. Zinn-Justin, Phys.Rev. D15 (1977) 1558- 
1564. 

[4] J. Zinn-Justin and U.D. Jentschura, Annals Phys. 313 (2004) 197-267. 

[5] J. Wosiek, Nucl.Phys. B644 (2002) 85-112. 

[6] M. Campostrini and J. Wosiek, Phys.Lett. B550 (2002) 121-127. 

[7] M. Trzetrzelewski and J. Wosiek, Acta Phys.Polon. B35 (2004) 1615-1624. 

[8] Piotr Korcyl, Acta Phys.Polon.Supp. 2 (2009) 623. 

[9] C.l M. Bender and T. T. Wu, Phys.Rev. 184 (1969) 1231-1260. 

[10] E. Brezin, G. Parisi, and J. Zinn-Justin, Phys.Rev. D16 (1977) 408-412. 

[11] C. Domb and M.S. Green, Phase Transitions and Critical Phenomena. 
3., Academic Press, London, 1974. 

[12] Z. Ambrozinski, Tunneling with Tamm-Dancoff method, Master's thesis, 
Jagiellonian Univ. Krakow, 2012 arXiv:1207.3898 [quant-ph]] 

[13] M. Unsal, [arXiv:1201.6426l [hep-th]]. 



9 



